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A two-dimensional, linear stability analysis of a thin nonplanar liquid sheet flow in vacuum is carried 
out A sheet flow created by a narrow slit of W and r attains a nonplanar cross section as a consequence 
of cylinders forming on the sheet edge under the influence of surface tension forces. The region where 
these edge cylinders join the sheet is one of high curvature, and this is found to be the location where 
Instability is most likely to occur. The sheet flow is found to be unstable, but with low growth rates for 
symmetric wave disturbances and high growth rates for antisymmetric disturbances. By combining the 
symmetric and antisymmetric disturbance modes, a wide range of stability characteristics is obtained. 
The product of unstable growth rate and flow time is proportional to the width-to- thickness ratio of the 
slit generating the sheet Three-dimensional effects can alter these results, particularly when the sheet 
length-to- width ratio is not much greater than unity. 


Nomenclature 

A = area; disturbance amplitude 

c = wave propagation velocity 

k ~ wave number 

L = length of sheet 

M = dimensionless parameter related to stability 
parameter 

R = radius of curvature 

r = x position at which end cylinder joins main sheet 

measured from the center of cylinder 
s = sheet-edge shape 

U - dimensionless x velocity 

w, v, w = velocity components in Cartesian coordinate 

system 

W = slit width 

We = Weber number 

X, K, Z = dimensionless parameters 

x, y, z - Cartesian coordinate system 

p = ratio of antisymmetric to symmetric disturbance 
amplitudes 

y = dimensionless parameter 

8 = stability parameter 

£ = dimensionless parameter related to stability 

parameter 

A = wavelength 

tj, $ = dimensionless Cartesian coordinate system 
p = density 

a = surface tension 

t = slit thickness 

<f> = velocity potential 

iff = angle between antisymmetric waves and z axis 

o) = wave frequency 
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Subscripts 

asy = antisymmetric 

c - end cylinder 

e = edge 

max = maximum value 

s = surface 

sy = symmetric 

x, y, z — partial derivatives with respect to, or components 
in, Cartesian coordinates 

£ tj, 6 = partial derivatives with respect to, or components 
in, dimensionless Cartesian coordinates 
0 = initial or steady-state value 

Introduction 

S TABILITY of incompressible, thin sheet flows has been 
of research interest for many years. Now there is renewed 
interest in sheet flows because of their possible application as 
low mass radiating surfaces. 1-4 Because of their low mass, 
near immunity to micrometeoroid damage, and simplicity, 
sheet flows are excellent candidates for a space radiator sys- 
tem. The objective of this study is to investigate the fluid dy- 
namic stability of sheet flows that are of interest for a space 
radiator system. A stability parameter is derived to relate the 
product of the disturbance growth rate and the flow time of 
the system with other parameters of the flow. 

Thin sheet flows are dominated by surface tension forces. 
As a result of surface tension at the edges of the sheet, a flow 
that begins with a dimension W perpendicular to the flow di- 
rection coalesces to a point at a distance Lin the flow direction. 
The resulting triangular sheet is ideal for a space radiator. A 
sketch of the geometry of a thin liquid sheet flowing through 
a narrow slit at z = 0 is shown in Fig. 1. Surface tension forces 
at the two edges of the sheet push the edges toward the z axis. 
As a result, as the flow moves in the z direction, the edge 
cross-sectional area A e grows. To satisfy mass continuity, the 
edges approach each other and Anally meet at the point 
In the following two sections, the Linear stability analysis 
will be developed. Following that, experimental results are 
compared to the theoretical results and an explanation for the 
discrepancy between the experimental and theoretical results 
is presented. 
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Linear Stability Analysis 

It is well known 3 that a small-diameter cylindrical liquid jet 
becomes unstable and breaks up into drops as a result of sur- 
face tension. This surface-tension-driven instability was first 
described by Rayleigh. 6 For the case of thin, planar liquid 
sheets, however, the linearized theory, which applies to small 
disturbances of the sheet boundaries, does not yield any un- 
stable solutions. Taylor 7 gives the solutions for two possible 
wave modes that may exist when a small disturbance is intro- 
duced into the flow. In the symmetric mode, displacements on 
opposite sides of the sheet are in opposite directions; in the 
antisymmetric mode, displacements on opposite sides of the 
sheet are in the same direction. Such waves are superimposed 
on the sheet cross section. Both wave modes arc shown in 
Fig. 2. 

If the liquid sheet is flowing into a gas rather than a vacuum, 
then the interaction of the gas and liquid will cause the sheet 
to break up. Squire 8 used a linear stability analysis to calculate 
the growth rate for thin sheets. Hagerty and Shea 9 considered 
finite thickness sheets. Later, Clark and Dombrowski,‘° using 
a second-order analysis, were able to calculate the breakup 
lengths of liquid sheets flowing into a gas. Recently, Rangel 
and Sirignano" completed a nonlinear stability analysis of this 
case. Crapper 12 obtained an exact solution for waves of arbi- 
trary amplitude on a fluid of unlimited depth when surface 
tension is the only driving force. He obtains a result for the 
maximum wave amplitude (amplitude/wavelength = 0.73). 

Stability of a Nonplanar Liquid Sheet 

Because of the sheet-edge shape, the sheet flows in this 
study are not planar. From Fig. 3 we can see that in the region 
where the constant thickness sheet joins the edge there is large 
curvature and, thus, large surface tension. Also, imperfections 
in the slits that arc used to form the sheet flows produce non- 
planar flow regions. It was expected that in these nonplanar 
regions unstable flow may occur. As a result, a linear stability 
analysis of a nonplanar sheet was carried out. 

Dispersion relations for symmetric and antisymmetric waves 
are derived in Appendix A. This analysis uses equations for 
perturbation of the boundary As and velocity potential A <t>> that 
are obtained from linearizing about the steady-state solutions 
s 0 and 4>q. Besides the usual neglecting of second-order terms 
in As and A<£, it was also assumed that the steady-state bound- 
ary velocities, uqj and v^, are constant. This was done so that 
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Fig. 3 Sample calculated cross section. 

the usual traveling wave solution is applicable. In this analysis, 
the flow has been assumed to be two dimensional because the 
sheet length is typically much greater than the width of the 
slit, and we consequently expect gradients in the z direction to 
be small. Cases where this assumption does not hold will be 
discussed in the section on experimental results. 

The dispersion relation for the symmetric wave is given by 
Eq. (A 18), and for the antisymmetric wave by Eq. (A 19). Here, 
o) ty is the complex frequency for symmetric waves, 10 ^ is 
the complex frequency for antisymmetric waves, and k is 
given by 

k = 2ttA (1) 

If either of the two roots for o) ty or <o My have positive imag- 
inary parts, then an unstable solution to Eq. (All) for the 
boundary s exists. As Eqs. (A 18) and (A 19) indicate, if Mo* = 
s 0 , x = 0, then both and < 0 ^ have only real roots, so that no 
unstable solutions exist. This is the result obtained by Taylor 7 
for the planar sheet. In this case, s 0 = r 0 /2 and the phase ve- 
locities for the two waves arc the following: 



For thick sheets (kr 0 /2 o°) the phase velocities are the same: 

c, y = c My - \ZlTTcripX (jkTo — ► oo) 


(3) 
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This agrees with the linear result given by Crapper 12 for t 0 — » 
o°. Crapper’ s 1 2 exact result for waves of arbitrary amplitude 
A is 


c % y = = V(2tto/pA)[1 + (tt 2 A 2 /4A 2 )] (To -4 *») (4) 

where A/A 0.73. Thus, for planar sheets, even for arbitrary 
amplitude, the phase velocity is closely approximated by the 
linear result 

For very thin planar sheets the phase velocities of the two 
waves differ significantly. For kr 0 / 2 — » 0, the antisymmetric 
wave phase velocity is the following: 


y = \/2aipr 0 = u t (jIctq -4 0) (5a) 

For the symmetric wave the phase velocity is the following: 

c #y = (27t/A)Vo’T 0 /2p = (ir7b/A)w, (t*t 0 ~4 0) (5b) 

It has been found 13 that antisymmetric waves propagate at the 
rate at which the sheet edge moves toward the z axis. This rate 
is called the edge velocity u t . 

For the symmetric wave the phase velocity is small for thin 
sheets. However, for the antisymmetric wave the phase veloc- 
ity is independent of the wavelength (no dispersion). This is 
similar to sound waves in a perfect gas or waves in an elastic 
string with tension 2a and mass per unit of length pT 0 , as was 
pointed out by Taylor. 7 Lee and Wang 14 used the string analogy 
in analyzing a thin annular liquid jet. 

Taylor 7 showed that the angle ip that the lines of constant 
phase make with the vertical is given by 

i/r = sin -1 (c/w 0 ) (6) 

where w 0 is the sheet velocity in the z direction and c is the 
wave phase velocity. Since c My = the lines of constant phase 
produced by antisymmetric waves must be parallel to the sheet 
edge, and this has been observed experimentally. 13 Since (ttt/ 
A) « 1, c, y » c M y , symmetric waves appear to be nearly 
vertical. Vertical waves, or striations, have been observed on 
the sheet as well. In general, the sheet is affected by both 
symmetric and antisymmetric disturbances. 

Now consider the nonplanar sheet where the steady-state 
boundary velocities Uq ^ and v 0 _, are finite. For the steady-state 
solution the boundary condition given in Appendix B applies 
[Eq. (B4)]. 

If we define 

X = &o, Y — o>'s 0 Iuq^ and y = a'/pu 2 0j s 0 (7) 
and neglect the slope in the y-z plane so that 

v 0j = (8) 

then we may rewrite Eqs. (A18) and (A19) as follows: 

Y% + isoJX tanh X)Y tr + yX\is 0 . x - tanh X) = 0 (9) 

Yly + is 0 J.X coth X)Y„ y + yX\is 0fX - coth X) = 0 (10) 

Before presenting the general solutions for Eqs. (9) and (10), 
consider the solutions for the very thin sheet flows, i.e., the 
limiting condition X — ► 0. From Eqs. (9) and (10) 

limT sy = 0 (11) 

lim = —jj 0 ..x (12) 

X-*0 


therefore 

lim = ^(vqj^o) (15) 

x-*o 

which implies 

lim 0 )^ y = lim Im{aCy} (14) 

x-*o x—c 

Thus, in the thin sheet limit there is no unstable symmetric 
solution. However, in the antisymmetric case the imaginary 
part of the frequency < 0 ^ is nonzero, and either a decaying- 
or a growing- wave solution is possible. Since As — exp[— iotf] 
~ exp[o>,/], Eq. (14) indicates that an unstable solution will 
result, if and only if, < 0 in the upper half of the sheet for 
ks 0 — ► 0. In the case of the sheet edge, this occurs in the region 
where the edge cylinder joins the sheet. 13 Experimentally, it 
has been found that holes develop in the sheet in this region. 

In the purely symmetric case, the sheet flow is stable in the 
limit of infini te disturbance wavelength, but is slightly unstable 
for finite wavelengths. For wavelengths on the order of those 
exhibited by typical experimental flows, the unstable growth 
rate is usually too small to produce sufficient growth for hole 
formation before the fluid flows into the point of coalescence. 
In the purely antisymmetric mode, on the other hand, the sheet 
is predicted to be highly unstable in the region where the end 
cylinder joins the sheet In the experimental flows studied, 
however, only occasional hole formation is usually observed, 
and sheets without holes can be formed. Hence, neither the 
antisymmetric nor the symmetric mode is sufficient alone to 
explain the observed sheet flow. Indeed, the nonlinear solutions 
of Clark and Dombrowski 10 for a liquid sheet in the atmo- 
sphere consist of both symmetric and antisymmetric modes. 
Mansour and Chigier 15 also observed both modes in their ex- 
periments. 

Stability of Combined Symmetric and Antisymmetric Modes 

If the sheet flow is affected by a disturbance consisting of 
both a symmetric and an antisymmetric mode, then the sym- 
metric mode has a damping effect that tends to reduce the 
unstable growth rate below that of the purely antisymmetric 
mode. A stability parameter 5 is defined as 

8 * t f fu = (L/woM = (W72iOo>, (15) 

where t f is the flow time (time for a fluid element to go from 
z = 0 to z = L) and h = 1 /o>, = l/Im{o/} is the growth rate of 
a disturbance. This parameter is used to measure the relative 
stability of the flow by indicating the amount of disturbance 
growth that can take place upstream of the point of coales- 
cence: a disturbance amplitude will increase by a factor of e* 
as the fluid travels from the slit to the point of coalescence. 
The stability parameter depends on three quantities: the sheet 
width-to- thickness ratio, the ratio of the sheet thickness to the 
disturbance wavelength, and the ratio of the antisymmetric am- 
plitude to the symmetric amplitude. 

To obtain we must derive the dispersion relation for com- 
bined symmetric and antisymmetric modes. From Appendix 
A, Eqs. (A16) and (A17) must be satisfied at the sheet-vacuum 
boundary. Solving Eqs. (A16) and (A17) simultaneously at the 
boundary (y = s 0 ) yields 

Y 2 + iso^pXY + yX\is^ “ P) = 0 (16) 

where X , K, and y are given by Eq. (7) and 

_ A^y + A ty tanh ks 0 
A tt y tanh ks 0 + A, y 


( 17 ) 
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It is straightforward to check that Eq. (16) reduces to Eq. 
(A 18) for the symmetric case when A„ y = 0 and to Eq. (A 19) 
for the antisymmetric case when A ty = 0. 

In a typical thin sheet flow, the wavelength is considered 
large compared to the thickness of the sheet. In the limit of 
long wavelength ( ks 0 « 1), the solution for o> to first order is 


If 5 is written in terms of the nondimensional variables of 
Eq. (20), incorporating co, = (uo a /s 0 )Yi and Y t = Im{T}, the 
following is obtained: 

S = ZXWir) (23) 


o/=^y-» 

*0 


-ikv 0 


( A gy ~f* A t yJCS() 
A. U ykSQ A ty 


(18) 


where Z, = Im{Z) and 


Z = (1 /r?){AlA c )Y 


(24) 


The growth rate of disturbance waves is the imaginary part 
of <o: 


Using Eq, (24) in Eq. (16), and expressing X and y in terms 
of dimensionless variables, yields 


<4 = lm{o>} = -> -kv 0j (±S-± (19) 

The solution is unstable if < 0, just as it is for the anti- 
symmetric mode alone, since here, a> t > 0. In the pure sym- 
metric limit, A uy — » 0, and so u) t — »■ —k 2 s 0 v 0t ,. Since ks 0 — » 0, 
a>i approaches zero and the stable symmetric solution (11) is 
recovered. Similarly, if A,y = 0, the disturbance is purely an- 
tisymmetric and a>i which agrees with Eq. (14), 

the growth rate of the antisymmetric mode alone. 

The parameter Y represents the growth rate of temporally 
varying disturbances. If this is large, i.e., the time constant of 
such growth is small compared with the flow time (6 > 1), 
then the sheet is said to be unstable because holes develop 
upstream of the point of coalescence. If the time constant is 
greater than, or on the order of, the flow time (8 < 1), then 
the sheet is said to be stable because no holes form upstream 
of the point of coalescence. Whether disturbances become 
large enough to produce holes depends in part on upstream 
fluctuations as well, so that a sheet with 8 somewhat greater 
than unity, but with only relatively small upstream perturba- 
tions, might still be free of holes. 

To compute the stability parameter of the sheet, it is nec- 
essary to solve for the shape and surface velocities of the sheet 
edge. This is done in detail in Chubb et al. 13 ; a summary of 
the relevant results is given here. In the course of the deriva- 
tions, the x component of velocity u is assumed not to depend 
on y. Furthermore, the edge-shape solution is based on a two- 
dimensional theory, but experiments have revealed that three- 
dimensional effects also influence the edge shape. 

To simplify the derivation, we introduce a number of di- 
mensionless parameters 


f- 


x 

7/2’ 


V 





where the subscript s indicates conditions at the surface bound- 
ary. Chubb et al. 13 showed that 

UJL& 6) = (1/7 j)[A(£ tyAjm (2D 


The function 17(f) was found to satisfy the following differ- 
ential equation: 




1 

2(1 + a) 


I*®’-* -£] 


(1 + iff* (22) 


Z 2 + 




(1 + *£)-”(«•% - j8) = 0 


(25) 


Solving Eq. (25) for Z, and taking the imaginary part, results 
in the stability parameter: 

(0 

where /3 is given by Eq. (17) and 

£= -ir(VnXiMc) 

M = [4tt 3 /(1 + v^mZ 2 ) 

The stability parameter depends on f, which is a function 
only of the dimensionless coordinate f. The largest value of 5 
occurs where f is a maximum. This maximum occurs very 
close to the point where the edge cylinder joins the plane sheet 
(f = r in Fig. 3). The maximum value was calculated based 
on the computed edge-shape solutions and is shown in Fig. 4. 
Note that Zm is a function only of the dimensionless edge 
cross-sectional area A c (independent of fluid properties), and 
that f is not a strong function of A c (except in a very small 
range of A c ) : More rigorously, fA c is approximately a linear 
function of A*, as demonstrated in Fig. 5. The results of Fig. 
6 indicate that = 1.4 is a good approximation over the 
range of A c of interest. 

Consider the solution to Eq. (26) for the sheet flows of ex- 
perimental interest. 1) rfk « 1 and 2) rfk « A^IA^ and 
rfk « A ly /A„y, so that /3 *** A^y/A^. (If A ny fA ty — 1, and 
condition 1 is satisfied, condition 2 occurs automatically.) 

Under these conditions, Eq. (26) can be expanded in a Tay- 
lor series. Retaining terms up to order (7/ A) 2 produces the fol- 
lowing result: 


(26) 


(27) 

(28) 


•4-="G)]G)(7)fe) <*> 


where a and 7 iR e are uniquely determined based upon A c . A 
is the dimensionless edge cylinder area integrated along the f 
axis and is a function of f. It is shown in Appendix B that A c 
is a linear function of z, and so A c is a linear function of 6. 
Equation (22) must be solved numerically; a representative 
edge-shape solution is shown in Fig. 3. 


Neglecting second-order terms in rfk yields the following 
for 8: 



7 « 1. fi - (30) 

A A,y 
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Fig. 5 Dependence of stability parameter on end cylinder. , 

true and , linear fit 


5 



$ 



5 



Using Can = 1.4, the maximum stability parameter is the fol- 
lowing: 




-« 1 * 
A 




(31) 


Recall that Eq. (30) and therefore, Eq. (31) are derived based 
on the assumption that t/A « both A m lA t y and A ly /A My . For 
the pure symmetric or antisymmetric case, it is necessary to 
go back to Eq. (26). Note that A„y = 0 for the purely symmetric 
case, so that in the long wavelength limit 


p 



7TT 

T 




(32) 


This applies only for < 0, since only in that case is there 
instability (co, > 0). For long wavelengths, 8^ — > 0, thereby 
confirming that the sheet is stable for long wavelengths. The 
value of Cty , defined in Eq. (32), has been found to be approx- 
imately a constant equal to 2.443. The result of Eq. (32) was 
compared with the results of the full equation (26) for various 
7 /A. As shown in Fig. 6, appreciable discrepancy emerges for 
7 / A = 2 X 10“ 2 , which is also where 8 > 1. 

To demonstrate the relative stability of sheet flow with re- 
spect to symmetric disturbances, consider as an example a 
sheet flow produced by a slit 3.44 cm X 75 fjm. Assume a 
pure symmetric disturbance with t/A % 4X 10~ 3 . Substituting 
into Eq. (32) gives 8 = 2.443 (W/t)(7 /A) 3/2 = 0.28: the flow time 
is well below the growth time constant and the sheet is ex- 
pected to be stable. 

In the antisymmetric limit, A^ = 0, so 





A 

TTTTf 


<5«y — ► 




( 33 ) 


6 



Fig* 6 Calculated stability parameter in pure symmetric mode. 

jfX - a) 1 x 10\ b) 4 x 10% c) 1 x 10% and d) 2 x 10 *. , 

exact solution [Eq. (26)] and — , approximate solution [Eq. (32)]. 

Since Wfr » 1 for a typical sheet, Eq. (33) confirms that a 
highly unstable sheet flow results if disturbances are purely 
antisymmetric. For example, in the case of the 3.44 cm X 75 
/xm slit, 8 = 130: the flow time is two orders of magnitude 
above the growth time constant. When Eq. (33) is compared 
with the result obtained by solving the full equation (26), very 
good agreement is found, and the error in the value of 8 is 
well below 1 % for the range of wavelengths of interest. 

A third special case exists if the symmetric and antisym- 
metric amplitudes are approximately the same. If A„y/A iy = 1, 
then p = 1 for all values of t/A, and Eq. (31) becomes 

W 7 T 

5^- I. 4-7 7 « 1, 0=1 (34) 

r A A 

For the previously considered slit example we obtain 8 = 2.6: 
the growth time is on the order of the flow time. From this we 
expect either stable sheet flow or only slightly unstable flow 
(occasional hole formation). Note that, if the vertical striations 
appearing on the sheet are the result of symmetric disturbance 
waves, then Eq. (34) is equivalent to 6 = 1.4n, where n is the 
number of striations on the sheet. Agreement of Eq. (34) with 
the result for 8 obtained by solving Eq. (26) is demonstrated 
in Fig. 7. Note that for the wavelengths of interest, agreement 
is quite good. The error in 6 , which arises as a result of the 
approximation Eq. (34), is about 10% when 7 /A = 2 X 10' 2 . 
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Fig. 7 Calculated stability parameter at p * 1. 


Comparison of Linear Stability Theory 
and Experimental Results 

Recent experiments have produced sheets that seem to defy 
the analysis of the previous section. In researching the edge- 
shape geometry and the stability of these liquid sheet radiators, 
two different types of facilities have been used. The first is a 
large vacuum facility that uses diffusion pump oil (Dow-Com- 
ing 705 silicone oil) as the working fluid. This facility, which 
was used in previous studies, 2 * 13 does not have view ports that 
allow a view of the entire sheet. As a result, it is not suitable 
for studying hole formation. However, hole formation resulting 
from instability predicted by linear theory has been observed 
in these long sheet flows. 

To view the entire sheet flow, a second smaller facility using 
water as the working fluid in air was constructed. This exper- 
iment used W = 20 cm. The t could be varied continuously 
with a resolution of 1 fim. The variable slit was attached to a 
large plenum that was connected to the city water supply 
through a control valve. A large plenum (plenum volume » 
sheet volume) is necessary to damp out pressure fluctuations 
produced by the water supply. If pressure fluctuations do occur 
upstream of the slit, holes will be produced in the sheet flow. 
Holes resulting from these fluctuations do not necessarily oc- 
cur at the sheet edges. The results from this experiment defy 
the analysis of the previous section. For large W7t, that theory 
predicts instability, and sheets form that are free of holes. 


Consider two sheets generated from the smaller experiment, 
with width-to- thickness ratios of 4000 and 4848 (length-to- 
width ratios of 1.6 and 1.2, respectively), both of which are 
completely devoid of holes. Theoretically, however, if we ap- 
ply Eq. (35) for the stability parameter pertaining to the com- 
bined symmetric/antisymmetric disturbance, with an estimate 
for t/A of 4 X 10“ 3 , we obtain 6 = 22.4 when WIt = 4000. At 
this value of 6, we expect the sheet to be highly unstable. In 
fact, even if we assume only symmetric disturbances and apply 
Eq. (32), we obtain 8 = 2.47 for W/r = 4000, implying an 
order of magnitude increase (e 247 = 12) in the amplitude of 
small disturbances by the time the coalescence point is 
reached. It is not reasonable to assume that only symmetric 
disturbances are experienced, especially since antisymmetric 
disturbances have indeed been observed. Instead, we must con- 
clude, at least in these cases, that our two-dimensional stability 
analysis is invalid. 

One possible explanation of the discrepancy between the 
two experiments is that the second facility generally produces 
sheets of significantly lower length-to- width ratios. For small 
L/W, three-dimensional effects are more important since the 
edge cylinders are growing more rapidly than for L/W » 1. 
Therefore, the two-dimensional stability theory results are 
more likely to agree with the large L/W results than the small 
L/W results. 

That three-dimensional effects become important for small 
L/W sheets is also suggested by measurements of the edge 
cylinder shapes for a sheet generated by the same facility. Fig- 
ure 8, which is taken from Chubb et al., 13 summarizes obser- 
vations about the shapes of the edge cylinders made by pho- 
tographing a sheet (length-to- width ratio of 2.62) from two 
different angles. Only the nearly circular shapes are predicted 
from a two-dimensional analysis (Fig. 3 and Chubb et al. 13 ). 
Although the length-to- width ratio of 2.62 seems to be suffi- 


y 




Initial, Nearly Circular 
Cross-Sectional Shape 


Flattened Out 
Cross-Sectional Shape 


Rebuilding or Peanut Like 
Cross-Sectional Shape 


Larger, Rebuilt, Nearly Circular 
Cross-Sectional Shape 


Fig. 8 Schematic diagram of edge cylinder variation in flow. 
t/A - a) I x 10" J , b) 4 x 10~\ c) 1 x 10“ J , and d) 2 x 10 _I . 
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cicntly larger than unity, it might be small enough that three- 
dimensional effects are important. 

If the three dimensionality of the problem were taken into 
account, the disturbance growth rates could be fundamentally 
different from those derived in the two-dimensional analysis. 
If the sheet shapes of Fig. 8 are the result of three-dimensional 
effects, then it is reasonable to conclude that the sheets having 
length-to- width ratios of 1.6 and 1.2 are also significantly af- 
fected by three-dimensional effects. On the other hand, length- 
to- width ratios for sheets generated by the larger facility are 
on the order of 10 or 12, an order of magnitude higher than 
unity, and we expect the analysis to be closer to the experi- 
ment. 


The sign in front of the last term on the left-hand side of Eq. 
(A2) is negative for the upper half of the sheet ( s > 0), and 
positive for the lower half of the sheet (s < 0). 

Now linearize Eqs. (A1-A3) by substituting the following 
expressions for <#>(*, y, r) and s(x, t): 

<t> = y) + A<f>(x, y, r) (A5) 

s = Sq(x) + A s(x, o (A6) 

where <fo> and s 0 are the steady-state solutions, and A <f> and As 
are small perturbations from the steady-state solutions. Sub- 
stituting Eqs. (A5) and (A6) in Eqs. (A1-A3) yields 


Conclusions 

In summary, the linear stability analysis for combined sym- 
metric and antisymmetric modes yields a result for 8 . The 
value of this parameter compares well with experimental re- 
sults using long (L/W ** 10) sheets in the larger facility. The 
symmetric mode alone predicts a 8 that is much too small, 
whereas the antisymmetric mode alone yields a 8 that is much 
too large. When both modes are combined in nearly equal 
amplitudes, the resulting value of 8 corresponds to a situation 
of occasional hole formation, as has been observed. 

Three-dimensional effects might explain, at least in part, the 
stability of sheets formed from high W/t slits, but having low 
L/W. This is offered only as a suggestion; more experiments 
would be necessary to determine whether decreasing L/W in- 
deed stabilizes the sheet, all else being held constant. A three- 
dimensional theoretical analysis of the edge cylinder shapes, a 
complex and tedious computational procedure to be sure, 
might provide additional insight. 

The results of this study, both theoretical and experimental, 
are nonetheless encouraging. Where the theory and experiment 
are not in agreement, sheets are fortunately more stable than 
the theory predicts. Thus, the stability parameter might be 
treated as an upper bound on the growth rate of disturbances 
on liquid sheets. Most important, it has been determined that 
stable sheets, free of holes and suitable for use in a space 
radiator, can be generated. 

Appendix A: Dispersion Relations for Small 
Amplitude Waves on Nonplanar Sheet 

We neglect the ^-direction gravitational force and any cur- 
vature of the sheet in the y-z plane. Therefore, assuming the 
liquid sheet is incompressible and iirotational the flow is gov- 
erned by 


d 2 <f> d 2 

+ ~dy 2 = 


0 


(Al) 


subject to the following conditions at the vacuum-sheet bound- 
ary. 

Bernoulli equation 


* 2 (A <t>) . d\A<f>) d 2 (f>o d 2 ^ 

+ -^“17 (A7) 

*4*) . a(A*> . S(A*,) <r (4 s)„ 

» + “- — + 'JoTo** 0 

(A8) 

*A s) d(A s) d(A 4>,) d(A 4>.) 

— r Uof — r So^c — : = 7 (A9) 

dt dx dx dy 

Here, Uo and v 0 are the steady-state velocity components, s 0j[ 
is the steady-state slope of the boundary in the x-y plane, and 
the subscript s denotes conditions at the boundary. In obtaining 
Eqs. (A8) and (A9), second-order terms in A <f> and As have 
been neglected and the steady-state solutions have been elim- 
inated. 

Now assume wavelike solutions for A <f> and As: 

A<*> = 4>(y)e Ktx -^ (A10) 

A* = &'*'—» (All) 

Substituting Eq. (A10) in Eq. (A7) yields 

^ - *** = 0 (A 12) 

This has the following solution: 

£ = y sinh ky + cosh ky (A13) 

where A„ y and A ty are constants. There arc two types of wave 
solutions to consider. In the symmetric wave solution (Fig. 3a), 
v = v, on the upper boundary and v = — v, on the lower bound- 
ary. For the antisymmetric case (Fig. 3b), v = v, on both the 
upper and lower boundaries. Assuming As « So, the following 
results are obtained. For the symmetric case 

A&y = (Acy cosh ky)e Kkx ~*^ (A14) 


+ ^ - Z h* 

dt 2 ( ‘ s) ~ p (1 + s])* 2 


= const 


(A2) 


Boundary velocity - fluid velocity 


ds ds dd> 

dr = Vx ^ Jt + = Ity 


(A3) 


where 


ds d 2 s 

and 


(A4) 


and for the antisymmetric case 

A^„y = (A uy sinh ky)e* kx ~"‘ ) (A15) 

Using Eqs. (A 10) and (All) in Eqs. (A8) and (A9) gives 

+ V aj 4>, + (o , '/p)k 2 S = 0 (A16) 

*>'£ = ks 0 _,$ + i$ y (A17) 

where $>, = d&dy, »' = <<> — ku,,^ and a' = cr( 1 + sl_ x )~ 3a . 
Substituting Eqs. (A14) and (A 15) and solving Eqs. (A 16) and 
(A 17) simultaneously at the boundary (y = s 0 ) y yields the fol- 
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lowing dispersion relations for the symmetric and antisym- 
metric cases. For the symmetric case 


where A, the sheet-edge cross-sectional area for 0 ^ x' < x 
(see Fig. 2), is given by 


(col y) 2 + &v 0 .,(tanh ks 0 )col y + k 3 (cr f /p)(is 0x — tanh ks 0 ) - 0 

(A 18) 


and for the antisymmetric case 

(coLy) 2 + ikv 0j (coth ks 0 )oolt y + k 3 {a' Ip){is 0x — coth ks 0 ) = 0 

(A 19) 


A(x,z) = 2 J s(x\z) dx' (B8) 

At the point x = r(z), where the end cylinder joins the sheet, 
we have A - A c and su = —(rf2)u„ where A c is the total sheet- 
edge cross-sectional area. Therefore, from Eq. (B7) 


dA c dA c 


(B9) 


In applying Eqs. (A 16) and (A 17) at the boundary, it was as- 
sumed that 5 = s 0 + Aj ** s 0 - Also, the steady-state boundary 
velocities and have been assumed constant, so that the 
assumed wave solutions with constant k and co are applicable. 

Appendix B: Continuity Equation Solution 

The continuity equation for steady-state incompressible flow 
is the following: 


du dv dw 

— H + — 

dx dy dz 


(Bl) 


Since the gravity force and z-direction surface tension force 
are being neglected, u„ w = h>o, and r = r 0 are constants. 
Therefore, from Eq. (B9) 

J fMi) 

s(x , z) dx = — z = T 0 zV2VKe (B 10) 

0 w 0 

w here W e is defined as (r/pwlr 0 , and u e is replaced by u t = 
V2cr/pr 0 . 13 
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where w, v, and w are, respectively, the jc-, y-, and z-direction 
velocities. Referring to Fig. 2, if Eq. (Bl) is integrated from 
0 ^ y ^ s(x y z), then Eq. (Bl) becomes the following: 


dsu 

a* 


dsw 

+ 

dz 


v. 


ds ds 

u, — - w, — = 0 
dx dz 


(B2) 


where u and w are average velocities defined as follows: 


w = - I u dv vv = - I w dy 
S Jo S Jo 


(B3) 


Since the flow is symmetric about the x axis, v (y = 0) = 0. 

At the sheet surface the fluid y-direction velocity v, must 
equal the surface speed ds/dt. Therefore, for steady-state con- 
ditions 


ds 


ds 


ds 


• = w — -b >v — = v, 
dr dx dz 


Using Eq. (B4) in Eq. (B2) yields the following: 
d(su) d(sw) 


dx 


dz 


= 0 


(B4) 


(B5) 


If the gravity force and surface tension force (curvature in the 
y-z plane) in the z direction are neglected, then w z « u„ 
where the subscripts denote partial differentiation. Therefore 


dsu _ ds 
dx W dz 


(B6) 


Beginning at x = 0, where s = 0 and u = 0, integrate with 
respect to x to obtain 


su = 


1 

2 


w 


dA 

dz 


(B7) 
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